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ABSTRACT 

We present deep Giant Metrewave Radio Telescope (GMRT) radio observations at 240, 330 and 610 MHz 
of the complex radio source at the center of the NGC 1407 galaxy group. Previous GMRT observations at 240 
MHz revealed faint, diffuse emission enclosing the central twin-jet radio galaxy. This has been interpreted 
as an indication of two possible radio outbursts occurring at different times. Both the inner double and dif- 
fuse component are detected in the new GMRT images at high levels of significance. Combining the GMRT 
observations with archival Very Large Array data at 1.4 and 4.9 GHz, we derive the total spectrum of both com- 
ponents. The inner double has a spectral index a = 0.7, typical for active, extended radio galaxies, whereas the 
spectrum of the large-scale emission is very steep, with a = 1.8 between 240 MHz and 1.4 GHz. The radiative 
age of the large-scale component is very long, ~ 300 Myr, compared to ~30 Myr estimated for the central 
double, confirming that the diffuse component was generated during a former cycle of activity of the central 
galaxy. The current activity have so far released an energy which is nearly one order of magnitude lower than 
that associated with the former outburst. The group X-ray emission in the Chandra and XMM-Newton images 
and extended radio emission show a similar swept-back morphology. We speculate that the two structures are 
both affected by the motion of the group core, perhaps due to the core sloshing in response to a recent encounter 
with the nearby elliptical galaxy NGC 1400. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (NGC 1407) — intergalactic 
medium — radio continuum: galaxies — X-rays: galaxies: clusters 



1. INTRODUCTION 

Quoting [Lara et aL] ( |1999| ), "radio sources are born, grow 
and finally . . . sleep ". The dormant phase in the evolutionary 
course of a radio galaxy and fate of the radio-emitting plasma 
after the cessation of the nuclear activity are among the most 
intriguing issues in extragalactic astronomy. The active stage 
of a powerful radio source, typically associated with an ellip- 
tical galaxy, can last several 10 7 up to few 10 8 years. During 
this time, the radio galaxy is likely fed by mass accretion onto 
the supermassive black hole of the host galaxy. Once this ac- 
cretion is interrupted, or becomes insufficient to sup port radio 
activity, the radio source enters a dying phase (e.g., Murgia et 
al.|2011| and references therein); the radio emission passively 
evolves and rapidly fades, even if expansion losses are negli- 
gible and the relativistic electrons are subject only to radiative 
losses. This is reflected into a pronounced steepening of the 
integrated radio spe ctrum whose slope a can reach ultra- steep 
values, a > 2 (e.g., Komissarov & Gubanov 1994), adopting 
the convention S v oc z/~ a for the synchrotron spectrum, where 
S v is the flux density at the frequency v. Eventually, the radio 
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emission disappears below the detection limit of present radio 
telescopes. 

During the fading stage of a radio source, the central nu- 
cleus may switch on again and produce new radio activity, 
thus leading to a restarted radio source. Evidence for episodic 
radio activities are reported in the literature for a growing 
number of radio galaxies (see, for instance, Saikia & Jamrozy 
2009 for a review) . Double-double radio galaxies (Schoen- 
make rs et al.||2 000) qualify as one of the most unmistakable 



examples of recurrent nuclear activity; here a new pair of in 
ner lobes are produced close to the nucleus before the pre 
vious ly generated, more distant ones have comple t ely faded 
(e.g., Subrahm anyan et aLl|1996| |Lara et al. 1999; |Saikia~& 
Jamro zy 2009). Another type of restarted radio sources are 



nesting radio galaxies (Jamroz y et al.| |2007), i.e., sources 
characterized by a faint, extended structure, within which a 
bright radio source, extended on a significa ntly smaller scale, 
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98% and 4C 29.30 ( [Jamrozy et al.|2007| ) 
cases, the large-scale emission and inner 
source appear clearly, and q uite sharply, separated in the spec- 
tral index distribution (e.g., Gizani & Leahy 2003), with much 
steeper spectral index values measured for the outer emission. 
This has been interpreted as evidence for different epochs of 
jet activity, with the more extended and steeper component 
being related to an earlier outburst of the central active galac- 
tic nucleus (AGN). 

In this paper, we focus on the complex radio source associ- 
ated with the elliptical galaxy NGC 1407 (see Table[j}, at the 
center of the poor group of galaxies Eridanus A at z=0.0059, 
with a one-dimens ional velocity dispersion a v =372 km s~ 



dBrough et al.|2006|) and .3-2 keV X-ray luminosity of 10 41 7 
Forbes et al.|2006| ). Based on Giant Metrewave Radio 
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TABLE 1 
General properties of NGC 1407 



RAj2ooo (h m s) 
DEC J2 ooo (° ' ") 
z 

D L (Mpc) 

angular scale (kpc/ /x ) 
morphology 
(mag) 



03 40 11.9 
-18 34 39 

0.0059 

25.0 

0.120 
E0 

-22.8 



Notes to Table [T]- RAj2ooo an d DECj2ooo are the coordinates of the optical 
galaxy from the NASA/IPAC Extragalactic Database (NED). The optical red- 
shift z is from Quintana et al. (1994). The galaxy morphological classification 



is from th e Third Reference Catalogue of Bright Galaxies (De Vaucouleurs etj 
20061. We assume a flat cosmology with Hq = 71 km s 1 Mpc L , fl\ ■ 



|1991) . The Cousins R-band absolute magnitude is from Trentham et al. 
anTO = 0.27. 



Telescope (GMRT) observations at 240 MHz and 610 MHz, 
Giacintucci et al. (2011, hereinafter Gil) suggested that this 
source may be a restarted radio galaxy. Here, we present a 
multifrequency radio study of NGC 1407 based on new, deep 
GMRT observations at 240 MHz, 330 MHz and 610 MHz, 
and multifrequency data from the Very Large Array (VIA) 
archive. Our study is complemented with the analysis of 
Chandra and XMM-Newton X-ray data. 

2. RADIO OBSERVATIONS 

In this section, we present the new GMRT observations of 
NGC 1407 at 240 MHz, 330 MHz and 610 MHz. Details on 
these observations are summarized in Table 2] which reports: 
observing date, frequency and total bandwidth in columns 1 , 
2 and 3; total time on source (column 4); full- width half max- 
imum (FWHM) and position angle (PA) of the full array (col- 
umn 5); rms level (la) at full resolution (column 6). 

The GMRT data were supplemented by all useful observa- 
tions in the VLA archive. The observing details are provided 
in Table [3] which shows the project code and array configu- 
ration in the first two columns; the other columns provide the 
same information as Table [2] 

2.1. GMRT observations 

NGC 1407 was observed with the GMRT at 240 MHz, 
330 MHz and 610 MHz (project 17_034; Table [2) in a full- 
synthesis run of approximately 9 hours (including calibration 
overheads) at each frequency. The observations at 240 MHz 
and 610 MHz were made using the dual-frequency mode, 
recording LL polarization at 240 MHz and RR polarization 
at 610 MHz. Both RR and LL were recorded at 330 MHz. 

The 330 MHz and 610 MHz data were collected using the 
upper and lower side bands simultaneously (USB and LSB, 
respectively), for a total observing bandwidth of 32 MHz. 
Only the USB was used at 240 MHz, with 8 MHz bandwidth. 
The default spectral-line observing mode was used, with 128 
channels for each band (64 channel at 240 MHz) and a spec- 
tral resolution of 125 kHz/channel. The datasets were cali- 
brated and reduced using the NRAO Astronomical Image Pro- 
cessing System package (AIPS). We refer to Giacintucci et al. 
(2008, 2011) for a description of the data reduction. 

Self-calibration was applied to reduce residual phase varia- 
tions and improve the quality of the final images. Due to the 
large field of view of the GMRT at all frequencies, we used 



the wide-field imaging technique at each step of the phase 
self-calibration process, to account for the non-planar nature 
of the sky. The final images were produced using the multi- 
scale CLEAN implemented in the AIPS task IMAGR, which 
results in better imaging of extended sources compa red to the 
traditional CLEAN (e.g., |Clarke & Ennslin |[20Q6j for a d e- 
tailed discussion, see Appendix A in |Greisen et al!~| p009)). 
We used delta functions as model components for the unre- 
solved features and circular Gaussians for the resolved ones, 
with increasing width to progressively highlight the extended 
emission during the clean. Beyond the image at full reso- 
lution (Table [2b, we produced images with lower resolution 
(down to ~ 45^, tapering the uv data by means of the param- 
eters ROBUST and UVTAPER in IMAGR. The final images 
were corrected for the GMRT primary beam response using 
the task PBCOR in AIPS. The rms noise level (la) achieved 
in the final images at full resolution is given in Table [2] 

We adopted the scale of |Baars et al.| ( |1977] ) for the flux den- 
sity calibration. Residual amplitude errors ar e within 8% at 
240 M Hz and 330 MHz, and 5% at 610 MHz ( [Chandra et al.| 
[2004] ). 

2.2. VLA archive data 

We extracted and reduced VLA archive data of NGC 1407 at 
327 MHz, 1.4 GHz and 4.9 GHz (see Table|3]for details). The 
observations were all pointed on NGC 1407, excepting those 
at 4.9 GHz in BnC configuration, where the phase center is 6' 
south-west of the galaxy. 

The observations at 327 MHz were obtained in multi- 
channel continuum mode with 16 channels and a total band- 
width of 6.3 MHz. The dataset was calibrated using 3C48 
as bandpass and amplitude calibrator. Due to problems with 
the source observed as phase calibrator (0340-213), 3C48 was 
also used for the phase calibration. 

All data were processed in AIPS and images were produced 
using the standard Fourier transform deconvolution method. 
Self-calibration was applied to reduce the effects of resid- 
ual phase errors in the data. Due to the large field of view 
at 327 MHz, the wide-field imaging technique was imple- 
mented to correct for distortions in the image caused by the 
non-coplanarity of the VLA. Correction for the primary beam 
attenuation was applied to the final images at all frequencies 
using the task PB COR in AIPS. All flux densities are on the 
Baar s et aL] ( |1977| ) scale and average residual amplitude errors 
are < 5%. 

3. THE RADIO IMAGES 

The GMRT 330 MHz image at full resolution (~12") is 
presented in Fig. HI, overlaid on the optical image from The 
Second Palomar Sky Survey (POSS 2). Faint, diffuse emis- 
sion (white contours) is detected around the central double 
source S3 (magenta contours), in good agreement with the 
240 MHz image in Gil (also shown in Fig.[8]). SI, S2 and 
S4 are discrete radio sources located, in projection, within the 
diffuse emission. SI and S2 have no optical counterpart on the 
POSS-2 image, while a very faint object seems to be associ- 
ated with S4, suggesting that all three sources are background 
radio galaxies. 

Hereinafter, we will refer to the two components of the 
NGC 1407 system as the inner double and large-scale emis- 
sion. 

3.1. The inner double 



Recurrent radio outbursts in NGC 1407 



3 



TABLE 2 

Details of the GMRT observations. 



Observation v Av t FWHM, p. a. rms 

date (MHz) (MHz) (min) ("x",°) (mJybearrT 1 ) 



Nov 19, 2009 


240 a 


8 


270 


16.1x10.9, 36 


0.25 


Jan 17, 2010 


330 


32 


330 


12.1x9.2, 42 


0.16 


Nov 19, 2009 


610 a 


32 


270 


8.2x4.4, 46 


0.05 



Notes to Table 2 - a: observed in simultaneous 240 MHz/610 MHz mode. 



TABLE 3 

Details of the VIA archive observations. 



Project 


Array 


Observation 


V 


Av 


t 


FWHM, p.a. 


rms 


code 




date 


(MHz) 


(MHz) 


(min) 


("x", ) 


(mJy beam -1 ) 


AS827 


B 


Apr 30, 2005 


322/329 


6.3 


130 


27.7x17.0,-16 


3 


AV151 


A 


July 3, 1987 


1446/1496 


12.5 


3 


2.2x1.3, 11 


0.1 


AR536 


CD 


May 28, 2004 


1385/1465 


50 


20 


76.2x20.8, 54 


0.08 


AW136 


BC 


Jun 23, 1985 


4835/4885 


50 


6.5 


5.1x2.0,71 


0.14 


AW112 


C 


Jun 2, 1984 


4835/4885 


50 


8.5 


15.9x4.4, 44 


0.07 



Fig. [I]} zooms on the central double source (S3). The region 
corresponds to the magenta box in Fig.[T^. The contours and 
gray-scale image are the emission detected at 610 MHz at full 
resolution (8" x A"). Thanks to the higher sensitivity, more 
extended emission associated with the double is visible in this 
new 610 MHz image, compared to the previous observations 
published in Gl 1 (also reported in Fig.[9j. 

The 2" -resolution image from the VLA-A data at 1.5 GHz 
(Fig. [lj, contours) shows the innermost structure of the 
source. The contours are overlaid on the VLA-BnC image 
at 4.9 GHz at resolution 5" x 2" . At ~1 kpc scale, the source 
appears a core-jet with position angle of ~30°, which is ap- 
proximately aligned with the axis of the double on larger scale 
(panel b). A faint blob of emission is visible to the east at 1.5 
GHz, apparently detached from the core-jet. The blob is not 
detected at 4.9 GHz. 

3.2. The large-scale emission 

The GMRT 330 MHz image, produced at the lower reso- 
lution of 45 // x45 // , is presented in Fig. [2j The new images 
at 240 MHz and 610 MHz are shown at similar resolution in 
Fig. 3 and Fig. 4a. The large-scale emission is detected at 
high levels of significance at all frequencies. Its total extent 
is ~80 kpc, consistent in all images. The brightest regions of 
this component are also visible in the VLA-CnD image at 1 .4 
GHz, restored with a circular beam of 60" (Fig. 4b). It is clear 
that the large-scale emission becomes progressively fainter 
with increasing frequency. The morphology also seems to 
change from the 1.4 GHz image, where the source shows a 
bipolar structure with east-west axis, to the 240 MHz image, 
where the extended emission becomes fairly amorphous and 
more extended in the north-south direction. This suggests that 
the radio spectrum of this component is very steep. 

4. RADIO SPECTRAL PROPERTIES, MAGNETIC 
FIELDS AND RADIATIVE AGE 



We studied the spectral properties - total integrated spec- 
trum and spectral index distribution - and derived the physical 
parameters of NGC 1407 using the Synage++ package (Mur- 
gia2001). 

1. We fitted the radio spectrum of the inner double and 
large-scale emission separately, and derived the injec- 
tion spectral index and break frequency z^reak for 
each component. 

2. We obtained a spectral index image of the system and 
extracted the distribution of a along the source structure 
as a function of the distance from the center. By fitting 
the observed spectral index trend, we derived ^break for 
the large-scale emission. 

3. We calculated the magnetic field B Qq assuming that the 
relativistic particle and magnetic field energy densities 
are in approximate energy equipartition. We assumed 
cylindrical geometry, a filling factor of unity and that 
there is equal energy in relativistic ions and electrons 
(k = 1). We also assumed that the magnetic field is un- 
ordered along the line of sight. We used the radio lu- 
minosity at 240 MHz, ai n j from the spectral modelling, 
and imposed a low-energy cutoff of 7 miri =100 in the en- 
ergy distribution of the radiative electrons, where 7 is 
the electron Lorentz f actor (e.g., |Brunetti et al.||1997" 



|Beck & Krause 2005). This allows us to take into ac 
count the contribution from relativistic electrons with 
an energy as low as ~50 MeV. 

Under a number of assumptions, the knowledge of the 
break frequency in the spectrum of a radio source al- 
lows us to estima te the time elapsed sin ce the source 
formation (e.g., Myers & Spangler 1985). We assumed 
that radiative (i.e., synchrotron and inverse Compton) 
losses dominate over expansion losses, and that the 
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Fig. 1 . — a: GMRT 330 MHz full resolution contours (white, magenta and yellow), overlaid on the optical POSS-2 image. The restoring beam is 12. 1" x 9.2" , 
p.a. 42°. The rms noise level is 0.2 mJy beam -1 . Contour levels are -0.6 (dashed yellow), 0.6, 1.2, 2.4 (white), 4.8, 9.6, 20 and 80 (magenta) mJy beam -1 . 
Labels indicate the discrete radio sources embedded in the extended structure, b: GMRT 610 MHz full resolution image (contours and gray scale) of the inner 
double (S3). The restoring beam is 8.2" x 4.4", p.a. 46°. The rms noise level is 0.05 mJy beam -1 . Contour levels are spaced by a factor of 2 starting from 0.2 
mJy beam -1 . Dashed contours correspond to -0.2 mJy beam -1 , c: VIA 1.5 GHz contours (FWHM= 2.2" x 1.3", lcr=0.1 mJy beam -1 ) overlaid on the VIA 
gray-scale image at 4.9 GHz (FWHM= 5.1" x 2.0", lcr=0.14 mJy beam -1 ). Contours scale by a factor of 2 starting from 0.3 mJy beam -1 . 



magnetic field strength is uniform across the source and 
remains constant over the source life time. We also as- 
sume that the relativis tic electron population is isotropic 
( Jaffe & Perol a|1973] ) and reacceleration processes can 
be neglected. We used the following equation 



trad = 1590 



d0.5 



[(1+zKrea^ 



-0.5 



to derive the radiative age £ ra d of both radio components, 
where t m & is expressed in Myr, z^break in GHz, and B eq 
and B C mb in fiG. B C mb = 3.2(1 +z) 2 is the equivalent 
magnetic field of the cosmic microwave background 



(CMB) radiation, i.e., the magnetic field strength with 
energy density equal to that of the CMB at the redshift 
z. 

Our analysis is presented in detail in the following subsec- 
tions, and results are summarized in Table [5] which provides 
the radio luminosity at 240 MHz, volume, observed spectral 
index, radiative model used for the spectral fitting, injection 
spectral index, break frequency, equipartition magnetic field, 
radiative age and equipartition pressure for each component. 

4.1. Integrated spectra 

We measured the flux density of the inner double and large- 
scale emission in similar regions at all frequencies, i.e., the 
extent of the double at 610 MHz (Fig. lb) and area covered 
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Fig. 2. — GMRT low-resolution image at 330 MHz. The image has been restored with a beam of 45.0" x 45. 0", p. a. 0°. The rms noise level is 0.5 mJy 
beam -1 . Contours start at +3 a and then scale by a factor of 2. The -3cr level is shown by the dashed contours. Labels indicate the discrete radio sources (Fig.Rl). 



by the diffuse emission at 240 MHz (Fig. 3). The flux den- 
sities are summarized in Table [4] along with the associated 
uncertainties (la) and angular resolution of the images used 
for the measurements. 

The VIA image at 327 MHz (not shown here) detects the 
inner double, but does not show significant diffuse emission 
around it, due to its lower sensitivity compared to the GMRT 
330 MHz image (Tab. 3). For this reason, Table [4] provides 
only the flux density at 327 MHz for the central double. 

We also inspected the VLSSj^] image at 74 MHz and found 
a weak source at the position of the inner double, with a 
peak flux density of ~400 mJy beam -1 . The source is not 
listed in the VLSS catalog, implying that its total flux den- 
sity at 74 MHz is less than the 5a local rms level (la=86 mJy 
beam -1 ). We therefore assume an upper limit of S74mh z =430 
mJy for the double. No evidence of significant emission over 
the large-scale structure is visible on the VLSS image. Given 
that the total area occupied by the diffuse component corre- 
sponds to -24 VLSS beams (each beam is 80" x 80"), we 
can place a 3a upper limit of 6 Jy at 74 MHz. 

In Fig. [5] we show the integrated radio spectra of the inner 
double (red points) and large-scale emission (black points). 



VLA Low-frequency Sky Survey (Cohen et"aT 2007 



The data points at 1.5 GHz and 4.9 GHz from A and BnC 
configurations (Tables 2 and [4]) are not shown in the plot; due 
to the high angular resolution and lack of short spacings, the 
images from these datasets detect only the innermost region 
of the double (Fig. lc). 

The overall spectrum of the inner double has a bs = 0.69 ± 
0.03 between 240 MHz and 4.9 GHz. A spectral steepen- 
ing seems to occur at high frequency, with a bs changing 
from 0.65 ± 0.05 below 1425 MHz to 0.75 ± 0.06 above. 
The large-scale emission has a much steeper spectrum, with 
a obs = 1.80 ±0.05 in the 240 MHz-1.4 GHz range. The up- 
per limit at 74 MHz suggests a low-frequency flattening of the 
spectrum, with a ^ s < 1.6 below 240 MHz. 

4.1.1. Spectral modelling 

We fitted the spectrum of th e inner double us ing a continu- 
ous injection (CI) model (e.g., Kardashev 1962), in which the 
source is continuously replenished by a constant flow of fresh 
relativistic electrons with a power-law energy distribution in 
a region of constant magnetic field. Under these assumptions, 
the radio spectrum has a standard shape, with a low-frequency 
spectral index representing the of the youngest electron 
population, and a high-frequency spectral index limited to be 
< ai n j +0.5 at frequencies above z^reak- 

The resulting best-fit for the inner double (red line in Fig. [5]) 
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Fig. 4. — a: GMRT low resolution contours at 610 MHz. The restoring beam is 41 .8" x 37.0. " , p.a. 17. The rms noise level is 0.15 mJy beam 1 . b: VLA-CnD 
low resolution contours at 1425 MHz. The image has been restored with a circular beam of 60 /; . The rms noise level is 0.08 mJy beam -1 . In both panels contours 
start at +3 a and then scale by a factor of 2. The -3 a level is shown as dashed contours. Labels indicate the discrete radio sources (Fig.[lJ?). 

yields a non-aged spectral index ^^=0.55^ ^ and ^break ~ netic field and radiative age of the central double are B eq ~7 
5.6 GHz. With these values, the estimated equipartition mag- jiG and t m & ~ 30 Myr (Table 5). 
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TABLE 4 

Flux densities of the NGC 1407 system 



NGC 1407 



component 



inner double (S3) . 



(MHz) 



S„ FWHM 
(mJy) ("x") 



large-scale emission 



74 
240 
327 
330 
610 
1400 
1425 
1.50 
4860 
4860 

74 
240 
330 
610 
1425 



<430 
270±22 
207±12 
230±18 
153±8 
87±4 
85±4 
37±2 
34±2 
23±1 

<6000 
945±76 
668±53 
194±10 
38±2 



80x80 (VLSS) 

15.4 x 12.5 (GMRT this work) 

28x17 (VLA this work) 

12.0x9.2 (FigJIJi) 

8.2x4.4 (FigJjjT) 

45x45 (NVSS) 

60x60 (Fig. 4b) 

2.2x1.3 (Fig.[lJ;) 

15.9x4.4 (VLA this work) 

5.1x2.0 (Fig.[lJ:) 



80x80 (VLSS) 
45.0x45.0 (Fig.fe 
45.0x45.0 (Fig.g 
41.8x37.0 (Fig, 
60.0x60.0 (Fig. 



47) 
A 7) 



1000 



100 



I 1 1 

large-scale structure 
a ob8 =1.80±0.05 

model ^ 0FF 
a inj = 0.7 (fixed) 

^break < ^30 MHz 

toirAt. < °- 6 




100 1000 
Frequency (MHz) 



Notes - The flux density of the large-scale emission has been obtained after 
subtraction of the contribution of the inner double (S3) and point sources SI, 
S2 and S4. 

The very steep spectral index of the large-scale component 
(a bs = 1.8) suggests that this emission is old. For th is reason, 
we modelled its spectrum usi ng a CIoff model (e.g., |Parma et| 
al.|2QQ7l|Murgia et al.|2011| ). This model assumes an initial 
phase of electron injection at a constant rate by the nuclear 
source (CI phase), followed by a switch-off of the nuclear ac- 
tivity. A dying phase then begins, with the radio emission 
rapidly fading subject to the energy losses of the relativistic 
electrons. It is assumed that the magnetic field strength is 
uniform within the source, energy losses are dominant with 
respect to other processes, and the pitch angle of the radiating 
electrons has an isotropic distribution with respect to the local 
direction of the magnetic field. We imposed an initial spectral 
index ai n j=0.7, similar to the observed a of the inner double 
(Tab. 5), and derived an upper limit of 230 MHz for the break 
frequency. The estimated magnetic field is B eq ~2.1 /iG and 
the total radiative lifetime is very long, i.e., t m $ > 320 Myr. 

The CIoff model also provides ^q ff/W. the dy ing to 
total source age ratio (for details see Par ma et al.|2007| ). We 
found that the active phase lasted for at least 128 Myr and 
the nuclear activity switched off not more than t ff ~ 192 Myr 
ago. 

4.2. Spectral index image 

We obtained images of the radio spectral index distribution 
by comparing sets of GMRT images at different frequencies, 
produced with the same cell size, uv range, restoring beam, 
and corrected for the primary beam attenuation. The images 
were aligned, the pixels with brightness below the 5a level 
were blanked, and finally the images were combined into a 
spectral index image. In Fig. [6^ we present our best spectral 
index image obtained by comparing the images at 330 MHz 
and 610 MHz, both restored with a circular beam of 45" tsl- 
dius and with noise levels similar to the images in Fig. 2 and 
Fig. 4a. We also show the distribution of the spectral index un- 
certainty based on the rms noise in both radio images (Fig. [6})) 
and overlay the 330 MHz contours on both images to provide 



Fig. 5. — Integrated spectra of the inner double (red circles) and diffuse 
component (black circles) using the flux density reported in Table [4] The 
solid lines are th e best-fit CI (red) and CIoff (black) models describee! in the 
text (Sect. |4~l~T) . 



a reference for the source morphology. We note that Fig. [6J) 
does not include the absolute flux scale uncertainty of 5% at 
610 MHz and 8% at 330 MHz, which results in a typical error 
on the spectral index of a a = 0.15. 

The inner double (S3), unresolved at this resolution, has 
an average spectral index of ~ 0.75 ± 0.02, consistent with 
a obs = 0.66 ± 0. 15 given by the flux densities measured at the 
two frequencies (Table [4]). Sources SI and S4 have both a ~ 
0.8, while the weak source S2 has a ~0.6. The large-scale 
component has spectral index values ranging between a ~ 
1.5 and a ~ 2.5, with uncertainties of ~ 0.06-0.08 in the 
brightest regions and up to ~0.25 at the edges. 

Fig. [7] shows the spectral index profiles along the eastern 
and western portions of the diffuse emission. We first mea- 
sured the flux density in circular regions on the individual im- 
ages at 240 MHz and 610 MHz (see inset), both restored with 
a circular beam of 45 "-radius, and then computed the cor- 
responding spectral index. The circles were placed following 
the brightest regions along the diffuse emission. The radius of 
each circle was set to 55" to ensure a sufficient signal-to-noise 
ratio and independent data points. In the plot we also show 
the data point corresponding to the inner doubl e S3 , obtained 
from the same images. As noticed in Section |4~T| there is a 
clear and abrupt separation in a between S3 and outer emis- 
sion, the latter being substantially steeper, with a = 1.5 ± 0.1 
to 2.2 ±0.1 along the western region and a = 1.8 ±0.1 to 
2.2 ± 0. 1 in the eastern one. 

We interpreted the spectral index trend in Fig. [7] in terms of 
radiative aging of the relativistic e lectrons by synchr otron and 
inverse Compton processes (e.g., |Parma et al.|2007| and ref- 
erences therein), assuming that radiative losses dominate over 
expansion losses and reacceleration. If the magnetic field and 
expansion velocity of the relativistic plasma are both constant, 
we can then estimate the minimum break frequency from the 
fit of the observed spectral index profiles. We used the re- 
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TABLE 5 

Spectral modelling results, equipartition magnetic fields and radiative ages 





^240 MHz 

(W Hz" 1 ) 


V 

(kpc 3 ) 


«obs 


spectral 
model 


C^inj 


Z^break 

(GHz) 




^ rad 

(Myr) 


toFFArad 


ta 
(Myr) 


*OFF 

(Myr) 


^radio 

(erg cm" 


" 3 ) 


inner double 


2.02 x 10 22 


140 


0.69 ±0.03 


CI 


55+ 014 


5.6 


7 


30 








2.8x10" 


-12 


large-scale 
emission 


7.07 x 10 22 


84 x 10 3 


1.80 ±0.05 


CIOFF a 
CIOFF b 


0.7 c 
0.7 c 


<0.23 
0.34 


2.1 


>320 
264 


<0.60 
0.65 


> 128 
92 


< 192 
172 


2.5x10" 


13 



Notes to Table 5- a: fit of the total integrated spectrum; b: fit of the spectral index trend; c: fixed. 
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Fig. 6. — Color scale image of the spectral index (a) and spectral index uncertainty (b) distributions between 330 MHz and 610 MHz. The image has been 
computed from images with a restoring beam of 45" x 45". Contours at 330 MHz and labels are as in Fig. 2. The spectral index uncertainty map is based on the 
rms noise in both radio images. The absolute flux calibration uncertainty is not included. 



lation ^break d~ 2 , where d is the distance from the core, 
which reflects the fact that the radiating electrons age as they 
travel away from the nucleus. The two profiles along the 
large-scale c ompon ent were fitted separately adopting a CIoff 
model (Sect. |4.1.1| ) with fixed to 0.7, i.e., the spectral in- 
dex of S3. The best fit, shown as blue lines in Fig. [7] gives 
a similar break frequency in the east and west portions of the 
source - ^ b reak = 330 ± 30 MHz and z^reak = 340 ± 30 MHz, 
respectively - and a ratio f ff Arad °f 0-6 and 0.5. Using the 
average between these values and a magnetic field strength of 
2.1 fiG (Table 5), we estimated a total age of ~ 264 Myr. The 
CI phase lasted approximately 92 Myr and the nuclear engine 
switched off ~172 Myr ago. 

The total age derived here is shorter than the lower limit of 
~320 Myr obtained from the fit of the total spectrum (Table 
5). A reason for such discrepancy may lie in the use of dif- 
ferent areas to extract the flux densities in the two cases; the 
region used to derive the total spectrum corresponds to the 
whole area occupied by the diffuse emission at 240 MHz, in- 
cluding the emission which is the most distant from the center 
and, thus, plausibly the oldest. The circular regions in Fig. [7] 
sample instead only a part of that region. 

Another consideration is that the spectral index between 
240 MHz and 610 MHz (a obs = 1.70 ±0.10) is slightly flat- 



ter than the average slope of the total spectrum, i.e., a b S = 
1.80 ±0.05 between 240 MHz and 1.4 GHz. This could result 
in a break frequency from the spectral trend fitting which is 
higher than that estimated from the total spectrum, thus yield- 
ing a shorter age (for the same magnetic field strength). 

Finally, it is important to bear in mind that both methods 
rely on the age determination from the break frequency, which 
has large uncertainties due to the n umb er of underlying as- 
sumptions (see discussion in Section [TT] ). 

5. X-RAY OBSERVATIONS AND DATA ANALYSIS 

NGC 1407 was observed during Chandra cycle 1, on 2000 
August 16, for -49 ks (ObsID 791), with the ACIS-S instru- 
ment operating in very faint telemetry mode. A summary 
of the Chandra mission and instrumentation can be found in 
Weisskop f et al.| ([2002 ). The data were reduced and analysed 
using CIAO 4.2 and CALDB 4.3 following techniques similar 
to those described in O' Sulli van et al.|p007| ) and the Chandra 
analysis thread^] The first half of the observation suffered 
from significant background flaring, and the final cleaned ex- 
posure time was 30.9 ks. 

Point sources were identified using the WAVDETECT task, 
with a detection threshold of lxlO -6 , chosen to ensure that 

10 http://asc.harvard.edu/ciao/threads/index.html 
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Fig. 7. — 240 MHz-610 MHz spectral index distribution as function of 
the distance from the center, computed using the 55 ;/ -radius circular boxes 
shown in the inset. The solid line is the best fit of the radiative model de- 
scribed in the text. 



the task detects <\ false source in the ACIS-S field, work- 
ing from a 0.3-7.0 keV image and exposure map. All point 
sources were excluded except the source corresponding to the 
AGN. Spectra were extracted using the specextract task. 
Spectral fitting was performed in xspec 12.7.0e. Abundances 
were measure d relative to the abundance ratios of lGrevesse &l 
Sauval| ( [l998l ). A galactic hydrogen column of 0.0542 x 10 : 



,22 



cm -2 and a redshift of 0.0059 were assumed in all fits. Spec- 
tra were grouped to 20 counts per bin, and counts at energies 
outside the range 0.5-7.0 keV were ignored during fitting. 

Background spectra were drawn from the standard set of 
CTI-corrected ACIS blank sky background events files in the 
Chandra CALDB. The exposure time of each background 
events file was altered to produce the same 9.5-12.0 keV count 
rate as that in the target observation. Very faint mode back- 
ground screening was applied to both source and background 
data sets. 

XMM has also observed NGC 1407 (Obsld 0404750101), 
with an exposure of ~66 ks. A detailed summa ry of the 
XMM mission and instrumentation can be found in lJansen~etl 
[aL] ( |2QQl( and references therein). The observation is centred 
close to NGC 1400, with NGC 1407 -9 ' off-axis. To pro- 
vide a comparison with Chandra image analysis, we reduced 
this data using SAS vlO following the methods described in 
IP' Sullivan eTaL] ( |20TTT ). The EPIC-MOS instruments were 
operated in full frame and the EPIC-pn in extended full frame 
mode, with the thin optical blocking filter. Periods including 
background flaring, when the total count rate deviated from 
the mean by more than 3a, were excluded, leaving useful ex- 
posures of ~38 ks (MOS) and ~26 ks (pn). Point sources 
were identified using EDETECT_CHAIN, and regions corre- 
sponding to the 85 per cent encircled energy radius of each 
source (except that at the peak of the diffuse emission) were 
excluded. 



6. X-RAY ANALYSIS 
6.1. The X-ray images 

We initially examined the structure of the hot intra-group 
medium (IGM) using exposure-corrected 0.3-2 keV Chan- 
dra and XMM-Newton images. Previous ROSAT observa- 
tions have shown NGC 1407 to be more X-ray luminous than 
NGC 140 0, and that there is a clu mp of emission between the 
two (e.g., [Mulchaey et al.| |2003 ). Figure [8] shows the XMM 
image, adaptively smoothed using the SAS asmooth task 
with a signal-to-noise ratio of 10. The image confirms the 
features observed in the ROSAT data. The 240 MHz emission 
(contours) extends well beyond the brightest part of the X-ray 
emission associated with NGC 1407. 

The clump of emission between NGC 1407 and NGC 1400 
consists of both point sources and diffuse emission. Although 
brighter to the east of NGC 1400, the diffuse emission extends 
across the galaxy and to the west as far as the edge of the 
field of view. The origin of this trail is unclear. NGC 1400 
is known to ha ve a velocity offset fr om the group mean by 
~1100 km s -1 ( Quintana et al.|l9 94), which is ~3 times the 
velocity dispersion of the group (cr v ~ 370 km s -1 ) and ~2 
times the sound speed v son ic (for a 1.2 keV gas, v son ic ~ 520 
km s -1 ). Therefore, the trail could be tidally or ram-pressure 
stripped material. However, there is no clear connection with 
NGC 1407. 

Figure^ shows the Chandra image, smoothed with a 10" 
Gaussian. The emission is asymmetric, being swept back 
from the northern side, forming two wings to the east and 
west, and more extended to the south. There is some indica- 
tion of an edge, or drop in surface brightness, on the northern 
side. To investigate this structure further, we modelled and 
subtracted the mean surface brightness distribution to reveal 
any residual structures. Two /3-models were used to represent 
the galaxy and group emission, and a flat model for the back- 
ground. Fitting was carried out in CIAO SHERPA, correcting 
for exposure using a monoenergetic 0.85 keV exposure map 
(the energy was chosen to match the peak photon energy of 
the data). After subtraction of the model, the residual map is 
smoothed to bring out any remaining structures. The result 
is shown in Figure [9j). An arc of positive residuals extends 
across the northern side of the galaxy, with strong negative 
residuals to its north and south. The arc is narrow in the cen- 
tre, and broader in the wings, particularly on the east side. 
This again suggests the presence of a surface brightness edge 
on the north side of the galaxy, with denser or cooler gas form- 
ing the bright arc, and the northern negative residual indicat- 
ing the sharp drop seen in Figure [9^. The negative residual 
in the galaxy core is likely caused oy the bright arc and cen- 
tral point sources biasing the inner part of the model to high 
values. We note that using an elliptical beta model does not 
result in any significant changes to the structures but fails to 
model the core and central point source properly. 

Figure [T0| shows the 0.3-2 keV exposure-corrected surface 
brightness profile extracted across the northern arc, using a 
125° sector, with the mean surface brightness profile across 
the remainder of the galaxy for comparison. The surface 
brightness drops further in the north than in the rest of the 
halo, and the profile first steepens (between ~5-10 kpc), then 
flattens outside ~10 kpc. There is no clear discontinuity in the 
profile, which would indicate the presence of a shock or cold 
front. However, the exposure may simply be too shallow to 
detect such a feature; relatively wide spatial bins are required 
to determine reliable surface brightness values. 
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Fig. 8.— XMM-Newton adaptively smoothed 0.3-2 keV image of the NGC 1407 group, with GMRT 240 MHz contours overlaid from Gl 1 (FWHM=48.5 // x 
31.9"). Contours start at 3a = 3 mJy beam -1 , and then scale by a factor of 2. NGC 1400, a group member, is marked, as well as the trail of X-ray emission 
associated with it. The black dashed arc marks the surface brightness edge in Fig. [9] 




Fig. 9. — a: 0.3-2 keV Chandra image of NGC 1407, Gaussian-smoothed with a 10" kernel, with overlaid the 610 MHz contours at the resolution of 
5.6" x 4.3" from Gil. Contours are spaced by a factor of 2 starting from 3cr=0.3 mJy beam -1 . The white arc marks the approximate position of the surface 
brightness edge, ~70" north of the galaxy core, b: 0.3-2.0 keV Chandra surface brightness residuals after fitting and subtracting two circular /3-models and a 
flat model for the background, binned by a factor 16 and smoothed with a Gaussian with a radius of 3 pixel (23.66"). Contours at 330 MHz, at the resolution of 
33" x 22", are overlaid starting at 1.5 mJy beam -1 and then scaling by a factor of 2. The white box indicates the area of panel a. In both panels the radio beam 
is indicated by the white dashed ellipse. 



6.2. Gas properties 



To determine the physical properties of the IGM, we ex- 
tracted a radial spectral profile centred on the galaxy. Since 
both the gas around NGC 1407 and the larger-scale group 
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Fig. 10. — Exposure-corrected 0.3-2 keV surface brightness in a 125° arc 
northward (black) and averaged over the remainder of the galaxy (grey). Er- 
ror bars indicate lcr uncertainty. 

emission show signs of disturbed, asymmetric structures, we 
chose to use simple circular annuli containing 1000-3000 net 
counts. The radio emission is observed throughout the core 
and over much of the Chandra field of view, therefore we do 
not exclude the corresponding region, since to do so would 
restrict the profile to the outer parts of the group. However, 
we note that if cavities are present coincident with the radio 
emission, this will affect the results of the spectral fits, both 
in terms of density (owing to the partial filling factor of the 
X-ray emitting plasma) and temperature (since the distribu- 
tion of temperatures along the line of sight will change). The 
spectra were fitted with a deprojected, absorbed APEC model. 
Abundances were tied between bins where necessary to sta- 
bilise the deprojection. The typical abundance was ~0.6 Z . 
Figure [TT] shows the resulting temperature and density pro- 
files, and the pressure profile derived from them. Pressure 
was calculated as P = nkT where n = 2n e . 

The temperature profile is i n reasonable agreement with 
previous spectral studies (e.g., Rasmussen & Ponman 2007; 
Diehl & Statler 2008), showing a cool core surrounded by 
hotter intra-group emission. Both density and pressure pro- 
files follow unexceptional approximate powerlaws. 

6.3. X-ray temperature map 

To examine the spatial variation of temperature in the gas, 
we pre pared a temperatur e map using the technique devel- 
oped by |David et al.| ( [2009] ) which takes advantage of the close 



correlation between the strength of lines in the Fe-L complex 
and gas temperature in ~1 keV plasma. Most of the emis- 
sion from such gas arises from the L-shell lines from Fe-XIX 
(Ne-like) to Fe-XXIV (He-like). For CCD resolution spectra, 
these lines are blended to form a single broad peak between 
approximately 0.7 and 1.2 keV. The centroid or mean photon 
energy of this peak increases with the temperature of the gas 
as the dominant ionisation state of Fe shifts from Fe XIX in 
0.5 keV gas to Fe XXIV in 1.2 keV gas, with Fe XVII and Fe 
XVIII providing the strongest line emission at the tempera- 
tures seen in the core of NGC 1407. The mean photon energy 
of the blended L-shell lines is independent of energy above 
-1.3 keV. 

We can thus estimate the temperature distribution of the gas 
by mapping the mean photon energy in the 0.7-1.2 keV band. 
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Fig. 11. — Radial profiles of temperature, density and pressure in 
NGC 1407. Density and pressure in the outermost bin may be overestimated 
since it includes emission from the outer parts of the IGM along the line of 
sight. 

We estimated the m ean gas properties from the spectral fits 
described in Section J&2| taking an abundance 0.6 Z as rep- 
resentative and simulated spectra for a range of temperatures, 
with redshift set to that of NGC 1407 and the hydrogen col- 
umn set to the Galactic value. These showed that for tempera- 
tures between ^0.5 and ~1.3 keV, the relationship is approxi- 
mately linear (JcT = -6.38+7.83<£>). Prior testing in other 
systems suggests that the relationship is relatively insensitive 
to variatio ns in abundance and hydrogen column ( |Q' Su llivan 
|et al.|2009j ). 

Using this technique, we produced the map shown in Fig- 
ure 12 The temperature in the galaxy core is 0.5-0.7 keV, 
rising to 1.2-1.3 keV north and south of the galaxy. We 
note that the highest temperatures in the map are somewhat 
lower than those in the radial profile (Fig.[TT]), since the map- 
ping technique becomes insensitive to temperature increases 
above ~1.3 keV. The cool core region is relatively circular. 
The black dashed arc marks the location of the X-ray sur- 
face brightness edge. The gas behind the edge seems to be 
cooler (~0.9-l keV) than the gas right in front (~1.2 keV), 
suggesting that the feature may be a cold front. However, the 
exposure is too shallow to allow temperature measurements 
on a scale comparable to the size of the feature. We there- 
fore cannot rule out the alternative possibilities that there is 
no strong discontinuity, or that a weak shock is present. A 
deeper Chandra pointing would be required to constrain the 
gas temperature behavior across the edge and determine the 
real nature of this feature. 

The cool (~1 keV) regions immediately east and west of 
the galaxy core correspond with the wings of the arc of excess 
surface brightness seen in Figure 9, showing that it is caused 
at least in part by the increased emissivity of this cooler ma- 
terial. This may indicate that material has been drawn out of 
the galaxy core and either heated or mixed with warmer IGM 
gas. 
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of the true source age. 

Our age derivation is also based on the assumption of uni- 
form and constant magnetic field across the source. A cal- 
culation of the effects of magnetic field evolution is beyond 
the purpo s e of the present paper, and we refer to [R udnick et 
al] ( [T994] ), |Jones et al.| ( [T999| ), jBlundell & Rawlings| (|2QQQ), 
and references therein, for detailed works on the implica- 
tions of magnetic field evolution for the source aging. In the 
case of very aged emission as that in NGC 1407, the rela- 
tivistic plasma may be partially mixed with the hot thermal 
gas. In this case, magnetized filaments can be produced in 
the radio volume due to plasma instabilities and the total syn- 
chrotron emission would result from the convolution of dif- 
ferent spectra produced by relativistic electrons emitting in 
regions with different magnetic field strength (e.g., Tregillis et 
al.|2004| ). The resulting convolution of the synchrotron kernel 



Fig. 12. — Chandra adaptively smoothed Fe-peak temperature map, with 
the low resolution 610 MHz contours overlaid (same as Fig.[4j?). The black 
dashed arc marks the surface brightness edge in Fig. [9] The radio beam is 
indicated by the white dashed ellipse. 



with magnetic field intensity and geometry yields a total spec- 
trum which is stretched and thus not straightfor wardly related 
to the spectrum of the emitting elect rons (e.g., |Eilek & A rendt 



7. DISCUSSION 



7.1. 



Restarted activity in NGC 1407 

The radio images presented in Section [3] show faint, dif- 
fuse emission, approximately 80 kpc across, which wholly 
encloses and dwarfs a small-scale (~ 8 kpc) double source at 
the center of NGC 1407. Based on former GMRT observa- 
tions, Gil suggested that the diffuse outer emission was pro- 
duced during an earlier cycle of activity of NGC 1407. The 
spectral analysis presented in this paper corroborates the mul- 
tiple outburst scenario. The large-scale emission is found to 
have an ultra-steep radio spectrum, with a = 1.82, as typically 
observed for highly evolved and fading radio sources (e.g., 
Komissarov & Gubanov|1994HParma et al.|2007[|Giacintucci| 



1996||Katz-Stone& Rudnic k[T997] ). In this case, standard ag- 
mg analyses, based on the position of the break frequency in 
the synchrotron spectrum, can give incorrect estimates of the 
age of the radio emitting electrons. 

7.2. Energy output of the radio outbursts 

We compared the energy output associated with syn- 
chrotron radiation from the two radio outbursts in NGC 1407. 
Using the radio luminosity at 240 MHz and in Table 5, 
we calculated the total radio luminosity over the frequency 
range 10 MHz - 100 GHz adopting the expression in |Q'Dea| 



& Owen| fL987] ) 



of Lto 



We found a bolometric radiative power 
1.6 x 10 42 erg s -1 for the inner double and L tot ~ 



3.2 x 10 42 erg s -1 for the large-scale emission. Based on our 
estimates of the radiative age of the two components (Table 
5), we found that the current outburst has released an energy 



of £tot,syn ~ 1.5 x 10 57 erg so far, while the energy of the for- 
et al.||2007[ |Murgia et a l. 2011). The spectral age analysis mer episode of activity, £ to t,syn ~ 2.6-3.2 x 10 58 erg, is one 

order of magnitude larger. 

It is well known that the synchrotron luminosity is an un- 
derestimate of the total energy output of a radio source, which 
is believed to be dominated by the mechanical work done 
by the radio jets on the external medium. The mechanical 
jet power of radio sources can be inferred, for instance, in 
those groups/clusters with depressions (cavities) in their X- 
ray surface brightness, interpreted as rising bubbles of rela- 



suggests that the radio plasma in such component is at least 
~300 Myr old. Furthermore, the radiative model that better 
describes the spectrum of this emission requires a switch-off 
of the nuclear engine ~ 170- 190 Myr ago, followed by a dy- 
ing phase. The inner double has instead a spectral index of 
a = 0.7 and a radiative age of ~ 30 Myr, consistently with 
being a currently active and relatively young radio source. 

The multi- scale radio morphology, combined with the dis- 
tinct spectral properties of its components, makes NGC 1407 
another remarkable example of a nesting radi o galaxy, other 



tivistic plasma inflated by the centra l AGN (e.g ., Birzan et al. 



nearb y examples being, for instance, 4C29 .30 ( Jamrozy et al. 
|2007] ) and Hercules A ( |Gizani et al.|2005| ). As double-double 
radio galaxies, these peculiar radio sources provide evidence 
for recurrent radio activity in elliptical galaxies. In the spe- 
cific case of NGC 1407, the small, young double source is 
currently fed by the central AGN, while the diffuse, steep- 
spectrum component is associated with relic plasma, which 
was injected during an earlier radio outburst of the AGN oc- 
curred at least ~300 Myr ago. 

It is important to note that our radiative ages have been es- 
timated neglecting adiabatic expansion. This could be a rea- 
sonable approximation for a relaxed, aged plasma in the dying 
phase such as that in the large-scale component. However, ex- 
pansion losses can be important for the inner double, as well 
as for the earlier, active phase of the large-scale component. 
Indeed, neglecting their effect may lead to an underestimate 



20041 |2008l ICavagnolo et aT1[20T0l [O^Sullivan et al.||2PTTT 
Here, it is estimated that the ratio of the mechanical to syn- 
chrotron luminosities is a few to few thousands for powerful 
radio sources, and up to several thousand for weaker sources 
(e.g., |Birzan et al.|2004| ). In those special systems with mul- 
tiple cavities, interpreted as signature of repeated AGN out- 
bursts, it is also possible to compare the energy outputs of the 
different episodes of activity. In some sys tems, such as Hy dra 
A ( |Wise et al.|2007) and NGC 5813 (Ra ndall et al.pOTT) , it 
is found that the past epoch of activity is the most energetic, 
suggesting that the mean jet power changed significantly over 
time or that the cu rrent outburst is stil l ongoing. In other sys- 
tems (e.g., A 262, Clarke et al. p009| Q, there is evidence for 
an opposite outburst trend, which may reflect an increase of 
the fueling of the AGN with time. 

The current X-ray images of NGC 1407 do not show evi- 
dence of X-ray cavities associated with the large-scale diffuse 
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radio structure. Dong et al. (2010) report a possible small 
(0.7x0.4 kpc) cavity in the group core, but its identification is 
dependent on the image processing employed and it appears 
to be uncorrelated with the small-scale active radio source. A 
direct measurement of the mechanical power of the radio jets 
in either period of activity is therefore not possible. If the 
mean jet mechanical power to radio power ratio has remained 
constant over the outburst history, NGC 1407 would then be 
another system in which the total energy output from the AGN 
is decreasing with time. On the other hand, the most recent 
episode of activity may be still ongoing (its radiative age is 
~30 Myr), and it is possible that the jet power may increase. 

7.3. Pressure comparison 

Assuming energy equipartition arguments, we can derive 
the non-thermal pressure in the radio source and compare it 
to the pressure of the surrounding X-ray gas. Under the as- 
sumptions listed at the beginning of Sect. |4| (see point 3), and 
adopting and B Qq in Table 5, we calculated a radio pres- 
sure in the large-scale component of P r adio,diffuse = 2.5 x 10~ 13 
erg cm -3 and P m dio S3 = 2.8 x 10~ 12 erg cm -3 for the inner dou- 
ble. 

Even though the X-ray pressure profile in Fig. 11 does 
not cover the whole extent of the radio emission, it is clear 
that the pressure of the X-ray gas is at least one order of 
magnitude higher than the radio pressure, with values rang- 
ing from ~ 10 -11 erg cm -3 at ~10 kpc from the center to 
~ 5 x 10~ 12 erg cm -3 at ~ 40 kpc. This is not unusual for 
cool-core systems, where a similar pressure discrepancy is 
often found (e.g., |Dunn et a l. 2005 ), suggesting a departure 
from equipartition conditions or an additional pressure sup- 
port in th e radio lobes, for insta nce, in the form of thermal 
gas (e.g., |Q' Sullivan et al.| |2010). Alternatively, an energet- 
ically dominant population of non-radiating relativistic par- 
ticles (protons) may contribute to the internal pressure. The 
achievement of pressure balance in NGC 1407 would then re- 
quire a ratio of energy in non-radiating particles to the energy 
in electrons k ~ 100- 1000. These values are in the range 
typically found for samples of centrally-located radio galax 



ies in groups and clusters ([Dunn & Fabian 2004 [ [Birzan et 
|aT1 [2008l |Croston et a^2008[ |Gitti et al.||2010j ) and suggest 
that that the lobes have been fed by heavy jets or have en- 
trained thermal material during t heir propagation throu gh the 
cluster/group atmosphere (e.g., Q' Sullivan et al. 2010). The 
diffuse and distorted morphology of the large-scale emission 
in NGC 1407, combined with the very long age inferred for 
this component, suggest that a mixing of the radio plasma and 
ambient thermal gas has already occurred at some level, fa- 
voring the entrainment scenario. 

7.4. Motion of the galaxy 

The morphology of NGC 1407 in the X-ray band (Fig. [8] 
and Figj9]) suggests that the galaxy is in motion. A surface 
brightness edge is visible in the north side of the galaxy with 
gas swept back to the east and west into two cool wings. There 
is some indication of a possible discontinuity in the surface 
brightness profile of the northern quadrant of the galaxy at this 
position (Fig. 10). If confirmed, this feature would suggest the 
presence of a shock or a cold front. The Fe-peak temperature 
map shows variations across the edge which seem consistent 
with the the hypothesis of a cold front, i.e., with the cooler 
gas behind the front (Fig. 12). However, the available data are 
not deep enough to constrain the temperature jump across the 



edge and thus the nature of this feature. 

Overall, NGC 1407 appears to be moving northward with 
its halo being stripped by the surrounding IGM. The location 
of the putative front is consistent with this scenario. If the 
edge is a cold front, it may indicate that the gas is sloshing 
in response to a recent interaction, such as a close passage, of 
NGC 1400. The X-ray trail visible in the XMM image (Fig. 8) 
could originate from the same interaction. 

The diffuse radio structure also shows a swept-back shape, 
but on a larger scale, similar to wide-angle tail (WAT) radio 
galaxies, although no jets are present. As discussed above, 
this emission is found to be the remnant of a former radio 
outburst of the central galaxy, consistent with the absence of 
jets. The X-ray wings seem to be inside the radio emission, 
and the radio contours appear compressed at the northern edge 

(e.g., Figjra. 

The similarities in the radio and X-ray morphologies sug- 
gest that the two structures are co- spatial and both affected by 
the motion of the galaxy. The X-ray wings may then represent 
a wake of cooler galactic material which has been stripped by 
the same ram-pressure that bent the diffuse radio emission. 

There are several ways in which we can obtain rough es- 
timates of the velocity of NGC 1407 relative to the IGM. If 
we assume, for simplicity, that the bending of the diffuse ra- 
dio emission is caused by ram pressure, we can use the Euler 
equation in the form 



F radio 



(e.g., O'Dea 1985 ), where r c is their curvature radius, r ra( jio 
is the radius or the radio tails, p ra di and v ra di are the density 
and velocity of the radio-emitting plasma, v ga i is the veloc- 
ity of the galaxy relative to the IGM, and /?igm is the den- 
sity of the IGM. From the images presented in Sect. [3] we 
estimate r rad i ~ 10 kpc and r c ~ 30 kpc. Based on the spa- 
tial extent of the diffuse radio emission (LLS~ 80 kpc) and 
age of the active phase tci ~ 90- 130 Myr (Table 5), we de- 
rive a first order estimate of the growth velocity of the radio 
source of v ra( jio = LLS/fci ~ 0.002c -0.003c, where c is the 
speed of light. Finally, from the Chandra data, we obtain 
Pigm ~ 0.0015 cm -3 w ithin the central 50 kp c and assume 
Pradio ~ 10~Vigm (e.g., Douglas s" et al.|[20lT ). This gives a 
very low velocity, only v ga i ~ 20 km s _1 . 

Alternatively, we can make a simple calculation assuming 
that the galaxy is moving north, that the radio outburst was 
located where we now see the southern boundary of the dif- 
fuse structure, and that it has moved at a constant velocity 
since that time. Neglecting projection effects and basing our 
estimate on the 240 MHz map, the galaxy would then have 
traveled 28-35 kpc over a period of ~300 Myr, suggesting a 
velocity of 90-115 km s -1 . 

These velocity estimates conflict to some extent with the X- 
ray morphology. If a cold front is present in NGC 1407, we 
would expect the galaxy to be moving at a significant fraction 
of the sound speed. For a temperature of 1 . 1 keV, as is found 
immediately north of the surface brightness edge, the sound 
speed IS V S onic ^ 

500 km s \ so a velocity of at least 250 km 
s -1 would be likely. 

However, there are large uncertainties on our estimates. If 
sloshing is taking place, the most likely cause is perturbation 
by a tidal encounter with NGC 1400. Given the 1100 km s" 1 
velocity offset between NGC 1400 and the group mean, it is 
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clear that sloshing motions would include a significant line of 
sight component. Since we see signs of a surface brightness 
edge in the X-ray and compressed contours in the radio, the 
velocity in the plane of the sky is greater than or comparable 
to that in the line of sight, but this still leaves an uncertainty 
of up to a factor y/2 in the true velocity of NGC 1407. The 
velocity could thus be as high as 160 km s -1 , or 0.3v ramc . 

Projection effects could have an impact on our estimate of 
the ram-pressure timescale. If the diffuse structure is in fact 
two old radio lobes, the axis of the jet which formed them 
may have been at an angle to the line of sight. In this case 
the lobes would be at a larger radius than we have assumed, 
and would experience lower external pressures. Mixing of 
entrained thermal plasma into the lobes would raise their ef- 
fective density. Both of these factors could increase the ram- 
pressure velocity estimate, though it is difficult to place limits 
on their impact. 

As discussed above, the X-ray images of NGC 1407 suggest 
that the group is not relaxed. A connection between bending 
of radio jets and sloshing induce d by minor mergers has been 
proposed in cluster cores (e.g.,|Clarke et al.|20 04; Ascasibar 
& Markevitch 2006 ; Mendygral et al. 2012 ), and it seems pos- 
sible that NGC 1407 is an example of this process at work in a 
galaxy group. In particular, we speculate that the group core is 
sloshing along the North-South axis, in response to a possible 
recent interaction with the nearby group galaxy NGC 1400. 
The gas motions induced by such sloshing may be then shap- 
ing both the X-ray and radio structures. There is a degree of 
tension between the velocities estimated from the radio and 
X-ray data. Resolving this issue would require deeper X-ray 
data, capable of detecting any cold front and providing a more 
reliable velocity measurement, and of determining the filling 
factor of the diffuse radio structure. 

8. SUMMARY 

We have examined the complex radio emission associated 
with the galaxy NGC 1407, at the center of the group Eridanus 
A. In a previous work, we have suggested that the galaxy ex- 
perienced two distinct radio outbursts. Thanks to new, deep 
GMRT observations at 240 MHz, 330 MHz and 610 MHz, 
combined with archival VIA data, we confirm the multiple- 
outburst scenario. NGC 1407 appears to be currently active in 
the radio band in the form of a small-scale (~ 8 kpc) double 
source associated with the galaxy. The analysis of the inte- 
grated radio spectrum indicates that the source has a spectral 
index of 0.7 and a radiative age of ~30 Myr, consistent with 
being an active and relatively young radio source. The double 
is wholly embedded in a faint, diffuse component, extending 
on a scale which is almost an order of magnitude larger (^80 
kpc). The large-scale emission is found to have an ultra- steep 
spectrum, with a = 1.8, as typically found for highly evolved 
and/or dying radio galaxies. The spectral age analysis shows 
that the radio plasma in such a component is at least ~300 
Myr old, supporting the idea that it was generated during a 
former radio outburst of NGC 1047. We estimated that the 
radio nucleus switched off ~ 170- 190 Myr ago, and after this 
time the source entered a dying phase. 

We compared the synchrotron energy output associated 
with the two radio outbursts and found that the energy re- 
leased during the duration of its current activity is nearly an 
order of magnitude lower than that from the former outburst. 
Numerous studies in the literature have shown that the syn- 
chrotron luminosity is an underestimate of the total energy 



output of a radio source, which is dominated by the mechani- 
cal work done by the jets on the external medium. If the mean 
jet power remained constant over the whole outburst history 
of NGC 1407, then NGC 1407 would be another example of 
a system in which the energy output is decreasing with time, 
other examples being Hydra A and NGC 5813. On the other 
hand, the most recent episode of activity may be still ongoing, 
and it is possible that the jet power may increase. 

We analyzed XMM-Newton and Chandra observations of 
the group. The diffuse emission surrounding NGC 1407 ap- 
pears swept back to the south, with wings extending to the 
east and west, and a possible surface brightness discontinu- 
ity on the northern side of the galaxy. Temperature mapping 
shows the emission inside the edge to be cooler than that to the 
north, but the data are shallow and the smoothing scale is too 
large to confirm whether a cold front or weak shock is present. 
The northern edge coincides with regions of compressed ra- 
dio contours, and the large-scale radio structure also has a 
swept-back morphology, indicating that the radio and X-ray 
structures are co-spatial and affected by the same forces. In 
general, the radio and X-ray morphology are suggestive of 
northward motion. We speculate that the group core may be 
sloshing on a north-south axis within the larger group poten- 
tial, having been tidally disturbed by an encounter with the 
nearby elliptical NGC 1400. Our estimated velocity in the 
plane of the sky is relatively low (~100 km s -1 ), but given 
the large uncertainties it seems possible that the true velocity 
could be a significant fraction of the sound speed, in which 
case sloshing motions could produce the radio and X-ray fea- 
tures we observe. 

Comparison of pressure estimates for the large-scale radio 
structure and the surrounding intra-group medium show the 
usual order-of-magnitude difference, the X-ray pressure es- 
timate being much higher than that found from the radio data. 
This may indicate a large population of non-radiating parti- 
cles within the radio structure, with a total energy ~ 100- 1000 
times that of the relativistic electron population. If, as seems 
likely, the structure was formed by a previous AGN outburst, 
these particles could be relativistic protons accelerated within 
the original radio jets, or a thermal plasma component which 
has been entrained during the propagation and expansion of 
the jets and lobes. The latter possibility provides an explana- 
tion for the continued presence of the radio structure in the 
group core, since buoyant forces, which would normally lift 
the lobes into the outskirts of the group halo, will be sup- 
pressed if the lobes contain a large fraction of thermal plasma. 
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